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Abstract 

The QCD radiative corrections to the leptonic decay rate of the Be meson are calcu- 
lated using the formalism of nonrelativistic QCD (NRQCD) to separate short-distance 
and long-distance effects. The Be decay constant is factored into a sum of NRQCD 
matrix elements each multiplied by a short-distance coefficient. The short-distance 
coefficient for the leading matrix element is calculated to order Os by matching a per- 
turbative calculation in full QCD with the corresponding perturbative calculation in 
NRQCD. This short-distance correction decreases the leptonic decay rate by approxi- 
mately 15%. 



The study of heavy quarkonium systems has played an important role in the devel- 
opment of quantum chromodynamics (QCD). Some of the earliest applications of perturba- 
tive QCD were calculations of the decay rates of charmonium [^. These calculations were 
based on the assumption that, in the nonrelativistic limit, the decay rate factors into a 
short-distance perturbative part associated with the annihilation of the heavy quark and an- 
tiquark and a long-distance part associated with the quarkonium wavefunction. This simple 
factorization assumption fails for P-wave states 0, which satisfy a more general factoriza- 
tion formula containing two nonperturbative factors [^. Calculations of the annihilation 
decay rates of heavy quarkonium have recently been placed on a solid theoretical founda- 
tion by Bodwin, Braaten, and Lepage Their approach is based on nonrelativistic QCD 
(NRQCD), an effective field theory that is equivalent to QCD to any given order in the 
relative velocity v of the heavy quark and antiquark p. Using NRQCD to separate the 
short-distance and long-distance effects, Bodwin, Braaten, and Lepage derived a general 
factorization formula for the inclusive annihilation decay rates of heavy quarkonium. The 
short-distance factors in the factorization formula can be calculated using perturbative QCD, 
and the long-distance factors are defined rigorously in terms of matrix elements of NRQCD 
that can be evaluated using lattice calculations. The general factorization formula applies 
equally well to S-waves, P-waves, and higher orbital-angular-momentum states, and it can 
be used to systematically incorporate relativistic corrections to the decay rates. 

Since the top quark decays too quickly to produce narrow resonances, the only heavy 
quark-antiquark bound states that remain to be discovered are the he mesons and their 
antiparticles. The possibility that he mesons may be discovered at existing accelerators has 
stimulated much recent work on the properties of these mesons and on their production 
cross sections at high energy colliders ||^. Once produced, a he meson will cascade down 
through lower energy he states via hadronic or electromagnetic transitions to the pseudoscalar 
ground state B^. which decays weakly. The discovery of the B^. meson will require a detailed 
understanding of its decay modes. In this paper, we compute the short-distance QCD 
radiative correction to the leptonic decay rate of the Be. We use the formalism of NRQCD 
to factor the amplitude for the decay into short-distance coefficients multiplied by NRQCD 
matrix elements. The short-distance coefficient for the leading matrix element is calculated 



1 



to next-to-leading order in as by matching a perturbative calculation in full QCD with the 
corresponding perturbative calculation in NRQCD. 

The leptonic decay of the proceeds through a virtual as in Figure 1. The 
W'^ couples to the Be through the axial- vector part of the charged weak current. All QCD 
effects, both perturbative and nonperturbative, enter into the decay rate through the decay 
constant fs^, defined by the matrix element 



where \Bc{P)) is the state consisting of a B^ with four-momentum P. It has the standard 
covariant normalization {Bc{P')\Bc{P)) = {2n)^2P'^6^{P' — P) and its phase has been chosen 
so that is real and positive. In terms of the decay constant fs^, the leptonic decay rate 
is 



where Vbc is the appropriate Kobayashi-Maskawa matrix element, Gi;' is the Fermi constant, 
Mb^ is the mass of the Be meson, and is the mass of the charged lepton. 

The formula ([l|) provides a nonperturbative definition of the decay constant Jb^, 
so that it can be calculated using lattice QCD simulations. One of the difficulties with 
such a calculation is that it requires a lattice with large volume and fine lattice spacing, 
since the strong interactions must be accurately simulated over many distance scales. The 
long-distance scales range from 1/Aqcd, the scale of nonperturbative effects associated with 
gluons and light quarks, to the scale l/{mcv) of the meson structure, where v is the typical 
relative velocity of the charm quark. The short-distance scales include the Compton wave- 
lengths l/rric and l/rrib of the heavy quark and antiquark. A more effective strategy for 
calculating /^^ is to separate short- distance effects from long-distance effects, to calculate 
the short distance effects analytically using perturbation theory in as, and to use lattice sim- 
ulations only for calculating the long-distance effects. Having already taken into account the 
short-distance effects, one can use a much coarser lattice which provides enormous savings 
in computer resources. 

An elegant way to separate short- distance and long-distance effects is to use NRQCD, 
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an effective field theory in which heavy quarks are described by a nonrelativistic Schroedinger 
field theory of 2-component Pauli spinors. The lagrangian is 

>Cnrqcd = /flight + ^1 {tDo + DV(2m,)) i;, + xl («^o - DV(2mfc)) Xb + ■ ■ ■ , (3) 

where £iight is the usual relativistic lagrangian for gluons and light quarks. The 2-component 
field ipc annihilates charm quarks, while Xb creates bottom antiquarks. The typical velocity 
V of the charm quark in the meson provides a small parameter that can be used as a non- 
perturbative expansion parameter. Lattice simulations using the terms given explicitly in 
@ can be used to calculate matrix elements for the Be with errors of order v"^. To obtain 
higher accuracy, additional terms, represented by ... in @, must be included. There are 8 
terms that must be added to decrease the errors to order v^. 

To express the decay constant Jb^ in terms of NRQCD matrix elements, we must 
express the axial-vector current Iry^'y^c in terms of NRQCD fields. Only the /i = component 
contributes to the matrix element (|1]) in the rest frame of the Be- This component of the 
current has an operator expansion in terms of NRQCD fields: 

^7°75C = Co{mb,mc) Xbi^c + C2{mb,mc) (D^fe)"^ ■ D^/'c + (4) 

where Co and C2 are short-distance coefficients that depend on the quark masses and m^. 
By dimensional analysis, the coefficient C2 is proportional to l/rriQ. The contribution to the 
matrix element {0\by^'~f^c\Bc) from the operator (DxbV ■ D-j/'c is suppressed by relative to 
the operator xli^c where v is the typical velocity of the charm quark in the B^. The ... in 
represent other operators whose contributions are suppressed by higher powers of f^. 

The short-distance coefficient Cq and C2 can be determined by matching perturba- 
tive calculations of the matrix elements in full QCD and NRQCD. A convenient choice for 
matching is the matrix element between the vacuum and the state \cb) consisting of a c and 
a 6 on their perturbative mass shells with nonrelativistic four-momenta p and p' in the center 
of momentum frame: p + p' = 0. The matching condition is 



(0|67°75c|c6) 



= Co {0\xl^c\cb) 

P-QCD 



+ C2 (0|(Dx6)T-D^,|c6) 

P-NRQCD 



+ (5) 

P-NRQCD 
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where P-QCD and P-NRQCD represent perturbative QCD and perturbative NRQCD, re- 
spectively. At leading order in a^, the matrix element on the left side of (||) is p)7°75Mc(p)- 
The Dirac spinors are 



«c(p) = \l ^"^J^" i pi c ] ' (6) 



2E. 



c \ Ec+iric 



where and t] are 2-component spinors and Eq = uiq + p^. Making a nonrelativistic 
expansion of the spinors to second order in p/mg, we find 

^.(-P)7°75«c(p) ^ vkc fi - I (!!^±±l!hL)\2 ^ \ 

At leading order in a^, the matrix elements on the right side of are rjl^c and p'^rjl^c- The 
short distance coefficients are therefore Cq = 1 and 

8 mj.gd 

where mred = ""^b^c/ (^6 + "^c) is the reduced mass. 

To determine the short distance coefficients to order as, we must calculate the matrix 
elements on both sides of to order a^. We will calculate the order-a^ correction only for the 
coefficient Co, since the contribution proportional to C2 is suppressed by v^. The coefficient 
Cq can be isolated by taking the limit p — 0, in which case the matrix element of (DxbV '^'^c 
vanishes. We first calculate the matrix element on the left side of @ to order a^. The relevant 
diagrams are shown in Fig. 2. The tree diagram in Fig. 2a gives a product of Dirac spinors 
Vb1o'j5Uc- We calculate the loop diagrams in Fig. 2b-2d in Feynman gauge, using dimensional 
regularization to regularize both infrared and ultraviolet divergences. Momentum integrals 
are analytically continued to D = 4 — 2e spacetime dimensions, requiring the introduction 
of a regularization scale /i. The effects of the quark self-energy diagrams in Fig. 2b and 2c 
is to multiply the matrix element by the renormalization constants y/Z^ and where 

= 1 + ^ -I^-9^-l(^°s47r-7) + ^log-^-l . (10) 
6t\ \ Aeuv ZeiR A 2 fi I 



The subcript UV or IR on the e's indicates whether the divergence is of ultraviolet or 
infrared origin. Combining the propagators using the Feynman parameter trick and using 
the fact that the external quarks b and c are on-shell, the vertex correction from Fig. 2d can 
be reduced to the tree-level diagram in Fig. 2a multiphed by the factor 
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After integrating over fc, it is convenient to change variables from the Feynman parameters 
to s = X + y and t = x/s. The s-integrals are trivial, but the t-integrals must be evaluated 
with care using the ie prescription in the denominator. Setting p-p' ^ m;,mc(l + f^/2), where 
V is the relative velocity of the b and c, and taking the limit v ^ 0, the integral reduces to 
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— 2 log h log 47r — 7 
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where mred = mbvud [j^h + ^c)- Multiplying the tree-level matrix element by the vertex 
factor 1 + A and by the renormalization constants \fZb and we obtain the final answer 
for the matrix element to order a,: 
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2 log h log 47r — 7 
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Note that the ultraviolet divergences have cancelled. The imaginary part of ([T3|) arises 
because it is possible for the quark and antiquark created by the current to scatter on-shell. 

We next compute the NRQCD matrix element on the right side of (|^) to order a,. 
The relevant Feynman diagrams are again shown in Fig. 2. The tree diagram in Fig. 2a gives 



a product of Pauli spinors r^^^c- It is convenient to calculate the loop diagrams in Coulomb 
gauge. In this gauge, the coupling of transverse gluons to heavy quarks is proportional to the 
heavy quark velocity, and therefore does not contribute in the limit v ^ 0. We therefore need 
only calculate the contribution from Coulomb exchange. The wavefunction renormalization 
factors associated with the diagrams in Fig. 2b and Fig. 2c are trivial: = = 1. The 
diagram in Fig. 2d reduces to a multiplicative correction to the tree-level matrix element 
from Fig. 2a: 



A 
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(27r)^ q2 + go - (P + q)V2m6 + - go - (p + q)V2me + le ' 

(14) 

where Eq = /2mQ. The infrared-divergent integral has been dimensionally regularized. 
After using contour integration to integrate over the energy go of the exchanged gluon, we 
find that (0) reduces to an integral over the gluon's three-momentum: 



^ _ 327ra^m,ed f ^'g 1 



3-2e n2 r,2 



+ 2p ■ q — ze 



(15) 



Evaluating the regularized integral in ([T5|) , we obtain 



P-NRQCD 
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TT^ ill / 
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+ log 47r — 7 



(16) 



where v is the relative velocity of the h and c: |p| = 'n^red'^^- Note that the infrared divergences 
and the factors of 1/f are identical to those in the matrix element (^) for full QCD. 

Comparing the matrix elements (|T3|) and ([16D , we can read off the short-distance 
coefficient Cq from the matching condition (^: 



Co 



TT 



■ log 2 

rrib + rric rric 



(17) 



We have chosen the scale of the running coupling constant. To the accuracy of this 

calculation, any scale of order or nic would be equally correct. Setting = 4.5 GeV 
and rric = 1-5 GeV, we find that asim^cd) ~ 0.34, and Cq ~ 0.85. 



6 



Our final result for the decay constant of the Be is 

ifsMs^ = Co {Olxli^cm + C2 (0|(Dx6)t-D^,|i?,). (18) 

The short-distance coefficient Cq is given to next-to- leading order in as in (|T^, while C2 is 
given to leading order in (^). The uncertainties consist of the perturbative errors in the short 
distance coefficients and an error of relative order f ^ from the neglect of matrix elements that 
are higher order in v"^. To achieve an error of order f^, the matrix element {OKDxbY ■'DipclBc) 
can be calculated by lattice simulations using those terms in the NRQCD lagrangian that 
are given explicitly in (^, but {0\xl'(pc\Bc) must be calculated using an improved action. If 
an error of order v"^ is sufficient accuracy, than the matrix element {0\(JDxbV ■'DipclBc) can be 
dropped, and {0\xl4'c\Bc) can be calculated using only those terms in the NRQCD lagrangian 
that are given explicitly in (|^). The parameters in this lagrangian are and the quark masses 
rrih and rric, all of which can be tuned so as to reproduce the spectroscopy of charmonium 
and bottomonium. These simulations should be very accurate for the intermediate case of 
the be system. In the absence of any lattice calculations of the be system, the matrix element 
{Olxli'clBc) can be estimated using wavefunctions at the origin from nonrelativistic potential 
models: 

mxli^c\Be)\' ^ 2M^ A|i?(0)p. (19) 

ZTT 

The factor of 2Mb^ takes into account the relativistic normalization of the state \Bc). 

Although it may not be relevant to the Be meson, it is interesting to consider the 
limit rrib ^ rric. In this limit, the perturbative correction in ([T7|) contains a large logarithm 
of mb/rric. Heavy Quark Effective Theory can be used to sum up the leading logarithms of 
rrib 0: 

\as{mh) J [ TT I rric \ j 

where ^/ is the number of light flavors, including the charm quark, and where A is an 
arbitrary factorization scale. 

We have calculated the short- distance QCD radiative correction to the leptonic decay 
rate of the Be meson. This calculation provides an illustration of the factorization methods 
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based on NRQCD that were developed in Ref. [Q. The decay constant Jb^ of the Be 
was factored into short-distance coefficients multiphed by NRQCD matrix elements. The 
radiative correction to the coefficient of the leading matrix element decreases the leptonic 
decay rate of the Be meson by about 15%. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40684. 
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Figure Captions 

1. Diagram for the annihilation of Be into lepton pairs via a virtual W~^. The shaded 
oval represents the wavefunction of the Be- 

2. Diagrams for the matrix elements (0167*^750106) in perturbative QCD and {Olxli'dci) 
in perturbative NRQCD. 
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